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ABSTRACT 

We report the detection of a ring of warm dust in the edge-on disk surrounding HD 32297 with the Gemini- 
N/MICHELLE mid-infrared imager. Our N' -band image shows elongated structure consistent with the orienta- 
tion of the scattered-light disk. The F v (ll.2fim) = 49.9 ±2.1 mJy flux is significantly above the 28.2 ± 0.6 mJy 
photosphere. Subtraction of the stellar point spread function reveals a bilobed structure with peaks 0."5-0"6 
from the star. The disk is detected out to the sensitivity limit at ~ 1", and the flux in each lobe is symmetric to 
within 10%. An analysis of the stellar component of the SED suggests a spectral type later than AO, in contrast 
to commonly cited literature values. We fit three-dimensional, single-size grain models of an optically thin dust 
ring to our image and the SED using a Markov chain Monte Carlo algorithm in a Bayesian framework. The 
best-fit effective grain sizes are submicron, suggesting the same dust population is responsible for the bulk of 
the scattered light. The inner boundary of the warm dust is located 0."5-0."7 (~ 65 AU) from the star, which 
is approximately cospatial with the outer boundary of the scattered-light asymmetry inward of 0."5. If the 
grains responsible for the asymmetry reside inward of the warm ring, they must have properties that differ from 
those that produce the bulk of the A^'-band emission. The addition of a separate component of larger, cooler 
grains that provide a portion of the 60 /im flux improves both the fidelity of the model fit and consistency with 
the slopes of the scattered-light brightness profiles. The models indicate an outer boundary of the small grain 
population beyond l."8 (~ 210 AU). Although the interpretation of the system is uncertain due to the unknown 
stellar age, previous indirect estimates (~ 30Myr) indicate the dust is composed of debris. The peak vertical 
optical depths in our models (~ 0.3-1 x 10~ 2 ) imply that grain-grain collisions likely play a significant role 
in dust dynamics and evolution. Submicron grains can survive radiation pressure blow-out if they are icy and 
porous. Similarly, the inferred warm temperatures (1 30-200 K) suggest that ice sublimation may play a role in 
truncating the inner disk. 

Subject headings: infrared: stars — circumstellar matter — planetary systems: protoplanetary disks — stars: 
formation — stars: individual (HD 32297) 



1. INTRODUCTION 

The attrition of the primitive remnants of solar system for- 
mation replenishes the dust in circumstellar debris disks. The 
scattered light and thermal emission from this dust provide 
a window into the physical processes governing the evolu- 
tion of solid material around normal stars, at a time after they 
have shed their primordial gas and dust envelopes and are 
transitioning to more mature, nearly dust-fr ee systems (e.g. 
iBackman & Parescdll993HMever et aT]|2007l). 

HD 32297 (d = 1 12 + \l pc: iPerrvman et al.lll997l) is a main- 
sequence star with a recently discovered circumstellar disk. 
ISchneider et all d2005l) first resolved the dust in near-IR scat- 
tered light with HST/NICMOS (F110W). They detected the 
near-edge-on inner disk out to 3. "3 (^400 AU), and found 
a brightness asymmetry inward of 0."5 and a break in the 
surface brightness profile a t 1"7. Optic al confirmation of 
the scattered-light disk by iKalasi (120051) followed, reveal- 
ing an extended, asymmetric outer disk extending to 15" 
(~ 168 AU) and sugg esting a blue-scattering R-[Fl 10W] dust 
color. iRedfieldl d2007l) detected the gaseous component of the 
disk, finding that HD 32297 exhibits the strongest Na I ab- 
sorption of any nearby main-sequence debris disk system. 

We report our discovery of spatially resolved emission 
arising from warm dust surrounding HD 32297 (^2|i, and 
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no te the independe n t disco very of resolved thermal emission 
by iMoerchen et aLl ((2007). We examine the morphology of 
the observed dust emission by removing the direct stellar con- 
tribution and determine that the residual structure suggests a 
model of thermal emission from an optically thin ring. We fit 
such a model to the data and explore the allowed parameter 
distribution (50. To conclude, we consider the implications 
of these findings in the context of physical processes shaping 
the disk's physical structure ($4| . 

2. OBSERVATIONS 

We observed HD 32297 and HD 20893 with MICHELLE 
on the Gemini-N telescope on the night of 2006 September 
19 (GN-2006B-C-12). HD 20893 serves as the photometric 
calibrator and point spread function (PSF) reference. We im- 
aged these stars in the N' filter (Ao = 11.2 /.im, AA = 2.4 /im) 
in a chop-nod sequence with a 15" throw. Observations were 
chopped and nodded at a position angle (PA) of 120°, cho- 
sen to be roughly perpendicular to the scattered light disk 
(PA = 47 .°6 ± 1°; G. Schneider, 2005, private communica- 
tion). Image quality was good (A^'-band resolution ~0"3), 
though terrestrial cirrus contributed to a varying background 
level which sometimes saturated the detector. The ABBA 
chop-nod sequences which contained saturated frames were 
discarded. Of the 188 s (23.5 s) of guided, on-source integra- 
tion time for HD 32297 (HD 20893), 94.0 s (23.5 s) was used 
in the final analysis. Each ABBA sequence was processed 
to remove both the sky and instrumental backgrounds via a 
double difference. The central, guided images resulting from 
each ABBA double difference were registered and stacked. 
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Photomet ric calibration wa s performed with HD 20893, 
which is a Coh en etail ([1999) standard with a zero-airmass 
mean A^'-band flux density of 4.29±0.14Jy at an isopho- 
tal wavelength of 1 1 .23 /Ltm. The mean airmass during 
the HD 32297 exposures was ~11% greater than during 
the HD 20893 exposures, and we assumed an extinction of 
0.172 mag AM -1 (corresponding to a mea n jV-band fAp 
1 1 ! I ! a m, AA = 5 fim] Mauna Kea extinction; Krisciunas et al.l 
1987) when correcting the measured HD 32297 fluxes. 

3. RESULTS & ANALYSIS 

The resulting images are shown in Figure Q] The image of 
HD 32297 exhibits extension consistent with the direction of 
the scattered light disk. The extended emission is detected 
out to the sensitivity limit at ~ 1". With a 1 ."4 radius aperture 
and an encircled energy correction derived from the image of 
HD 20893, we measured the total flux of the star and disk to 
be 49.9 ± 2.1 mJy. The uncertainty includes the contribution 
from both the background noise and the zero point uncertainty 
derived from HD 20893. As the filter is relatively narrow, a 
color-correction to this monochromatic flux density was not 
considered necessary. We place the measurement of total flux 
in relation to the known SED in Figure[2] 

3.1. Stellar Properties 

In order to characterize the dust emission, we must (1) 
estimate the fraction of the observed N' flux attributable to 
the star, and (2) estimate the stellar luminosity, which af- 
fects the temperature balance of the grains. To these ends, 
we characterize the star by modeling the optical and near- 
IR portion of the S EP, g i ven by Tycho-2 and 2MASS pho- 
tometry dHH^etal] HoS [StotskieitinilOOi). For model 
comparison, we synthesized photo metry with the stellar at - 
mospheres of the NextGen grid (lHauschildt et all fl999h . 
The scaling for the resulting photospheric fluxes, parameter- 
ized by £ = (R*/d) 2 , was least-squares fit to match the syn- 
thetic photometry to the optical/near-IR data. The choice of 
model atmosphere parameters is unclear at first glance, since 
the spectral type of HD 32 297 is inconsistent in the literature. 
It was classified as AO by ICannon & Pickering! in the Henry 
Dra per Catalog, while it specified as A5 in the AGK3 Cata- 
log dHeckmannl 1 1 9751) . Assuming no interstellar extinction, 
logg = 4.5, and [Fe/H] = 0, we found that a T eS = 9600 K 
spectrum (appropriate for A0V) is clearly rejected by the pho- 
tometry (xl = 200), while a 7600 K NextGen model fit the 
data best (xl = 2.5). The resulting best-fit model is shown in 
Figure |2] 

Can interstellar reddening bias our inferred T e ff? At 
^112pc, the star m ay be ou t side th e relatively dust-free Lo- 
cal Bubble, which Redfield (2007 ) notes is expected t o ex- 
tend to ^90 pc in this direction (Lallement etd]l2003h . As 
an empirical test, we fit for the extinction (par ameterized 
by A v ) using the A\/A v law of iFitzpatrickl (120041) and as- 
suming Ry = 3.1. The best-fit model, fixing T e ff = 9600 K, 
gives Ay = 0.72 ± 0.02 mag (formal error) and is rejected with 
xt = 8. 5. 3 When adjusting T & s, the 7600 K model again min- 
imizes x 2 7 with a best-fit Ay = 0.026 ±0.024 mag at this tem- 
perature. Thus, we empirically determine that the observa- 
tions are consistent with zero interstellar extinction. 

As an alternative approach to inferring the amount of in- 
terstellar dust toward the star, we consider the measurement 

3 With the addition of A v as a parameter, the number of degrees of freedom 
v has decreased relative to previous fits. 



of the interstellar component of Na I absorption by iRedfieldl 
(20071). He measures a column density of A^ai ~ 10 cm . 
There is considerable scatter in the correlation between A^jai 
and E(B-V), and assuming Ry = 3.1, values of Ay = 0.03 to 
0.3 mag are reasonable (Fig. 4d of lHobbsl 19741) . Since this is 
consistent with a small amount of extinction, we do not alter 
our conclusions. 

We estimate the flux scaling £ = (7.22 ±0.14) x 10" 20 from 
the best-fit zero-extinction model (after scaling the formal 
uncertainty by yxf). Using this quantity, we calculate the 
monochromatic flux density /v(11.23/im) = 28.2 ± 0.6 mJy. 
The uncertainties in our estimates of the stellar radius and lu- 
minosity are dominated by distance errors, as 

U = 4TtCd 2 aT^. (1) 

Taking the best estimates £ = 7.22 x 10~ 20 and d = 1 12 pc, we 
find R* = 1. 34 R@ and L» = 5.4L Q . Together, these suggest 
the literature value of A5 is closer to the "true" classification, 
though the cool temperature we estimate suggests even later 
A subclasses. As noted by Schneider et al. (2005), the star 
lies near the bottom of the A star main sequence in a color- 
magnitude diagram, a property in c ommon with other young 
debris disk systems (lJura et al.lll998l) . Figure 2 of iMoor et al.l 
(2006) shows that HD 32297 is relatively under-luminous and 
is positioned toward the red end of the local A star sequence, 
consistent with our results. A modern high-resolution spec- 
trum can settle the debate over the correct stellar classifica- 
tion. 

Finally, we note the ultraviolet measurements from the TD- 
1 satellite (Thompso n et all"l978l) shown in Figure|2] We did 
not use these data in our photometric fit. The F1965 and 
F2356 bands (A = 196.5 and 235.6 nm, respectively) show 
excess above the model photosphere. These may be due to 
line emission — such as Si III] and C III] in the former band, 
and Fe II in the latter — similar to what is seen in T Tauri 
stars (e.g. lLamzinl l 20 00a; Val enti et al.ll2003l) . In such stars, 
the intercombinational lines are lik ely associated with accre- 
tion rather than the chromosphere (Lamzin 2000b). The UV 
excesses i n HD 32297, a long with the recent evidence for a 
gas disk (Redfield 2007), lend indirect support to the con- 
tention that this is a young stellar system. 

3.2. PSF Subtraction 

Using the range of £ estimated in the previous section, we 
scaled the image of HD 20893 and subtracted it from the that 
of HD 32297. After PSF subtraction, we find 21 .6 ± 2.2 mJy 
in non-color-corrected residual flux (Fig. |3). The residu- 
als suggest the observed image (Fig. QJ) is a composite of 
a bilobed structure and an unresolved central source, whose 
flux can be fully accounted for by the star. The lobes peak 
at offsets of 0"5-0."6 from the star, corresponding to a ra- 
dius of ~60 AU. This suggests an inner edge to the population 
of grains with optically thin A^'-band emission — both hotter 
dust closer to the star and the case of optically thick material 
would fill in the emission at smaller projected separations. 

We quantified the degree of asymmetry between the disk 
lobes with aperture photometry of the PSF-subtracted disk 
image (Fig. |3J;). We used 0."6 square apertures placed on 
each lobe, with centers 0"6 from the stellar centroid. Ne- 
glecting an aperture correction and calibration uncertainty, we 
measured a NE-SW flux difference of 0.17 ±0.36mJy. As a 
fraction of the total flux in these apertures, the difference is 
(2.0±4.3) x 10~ 2 . We conclude that the disk emission is con- 
sistent with symmetry. 
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FIG. 1. — Panel (a) displays the result of our Gemini-N/MICHELLE imaging of the HD 32297 debris disk (N' band), showing elongated emission around 
the star. Contours are spaced according to the 1-ct background noise level. The elongation PA is consistent with the scattered light disk. The chop direction 
is indicated by the grey line on the compass. Panel (b) shows the image of the reference star HD 20893. The contour spacing is selected to match those in (a) 
relative to the stellar flux. 
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FIG. 2.— The SEP of t he HD 32297 system. IRA S photometry is from IMoor etafl J2006I) (see also [S ilverstone 2000). Near-IR data is from 
2MASS <Skrutskie etH]|2003). opt ical data from Tycho-2 IHog et alj|2000h . and UV data from TD-1 IThompson et al.ll 19781) . A 7600 K NextGen model 
photosphere (Hauschildt et al. 1999) with zero foregrou nd e xtinction was fit to the optical-NIR photometry. The range of allowed dust emission from Model I is 
shown in blue, while that of Model II is shown in red ( j|3.31 . Model I fits for the image and all SED points, while Model II treats the 60 fim measurement as an 
upper limit, effectively allowing for a separate population of cool dust. 
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FIG. 3. — Images of the disk with the stellar PSF subtracted for various 
values of the stellar flux estimate. With F® = 28. 2 mJ y and a = 1.1 mjy 
(which includes both error in our flux estimate, 41 3 - 1 1 and N' calibration 
uncertainty, ^2), panels (a)-(e) depict subtractions of a central source with 
F v = {F°-2a,F^-a,F$,F° + a,F° + 2a}, respectively. As in Fig.Q] con- 
tours are spaced according to the 1-cr background noise level. The bilobed 
residuals suggest that the system may be modeled as a central star with an 
optically thin ring of warm material. 

3.3. Modeling 

3.3.1. Model Construction 

We adopt the model for th ermal emission from o ptically 
thin dust rings developed by Backman et al. (1992), which 
assumes particles of effective size Ao that radiate as modi- 
fied blackbodies with emission efficiencies e„. We model the 
three-dimensional emission of the dust j u (x) by evaluating 



Mx)= r ° 
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jJx) = aJx)e^B l ,[T(r)]. 



(3) 

(4) 
(5) 

Here, a„ is the absorption coefficient, and the function f(z, h) 
describes the vertical distribution of dust and has scale height 
h. The fiducial values for the vertical optical depth to absorp- 
tion ro and scale height ho are set at the inner edge of the dust 
annulus, G3o. We use cylindrical coordinates for positions x, 
with r 2 = C3 2 + z 2 , and assume azimuthal symmetry. 

Some assumptions of grain properties are implicit in Equa- 
tions[2]-[5] The disk is assumed to be optically thin at all radii. 
The parti cles are efficien t abso rbers, but inefficient emitters 
(Eq. [3}. Backman et al. (1992) discuss the relationship be- 
tween Ao and the particle absorption efficiency and size dis- 
tribution in their Appendix D. They argue that for "dirty ice" 
dielectric consta nts and dn(a) oc a~ 7 / 2 da collisional-cascade 
size distribution ( Dohnanvil l 1 9691) . the minimum grain size is 
an uncertain a m j n ~ Aq/6. With these assumptions, the model 
contains only minimal information regarding grain character- 
istics. While we ignore the effects of grain composition with 
this approach, the absence of evidence for grain mineralogy 
(e.g. mid-IR spectroscopy) renders consideration of various 
compositions premature. Future work using more advanced 
models with explicit dependence on grain sizes and mineral- 
ogy may produce significant differences from our results. 

For simplicity, we assume the vertical density distribution 
/ is Gaussian, with width a = ho = 10(ii/112pc) AU fixed 
at all radii (rj = 0). Bec ause we do not measure a significant 
brightness asymmetry Q3.2b . we assume the disk is spatially 
centered on the star. The distance d is a 'nuisance' parameter 



not directly related to the disk architecture, and in the inter- 
est of independence we define 9q = XS^/d (likewise for 9\). 
The remaining free parameters are 8q, 9\, Ao, vertical optical 
depth to absorption to = t* s (GJo), (subpixel) stellar centroid 
x*, stellar flux (parameterized by £; jj3.ll >. disk inclination 
i, and PA. The stellar luminosity is given by Eq. [T] There 
are 1 1 degrees of freedom within the set of model parameters 

6 = {d,x,,PA,z',£,7,0o,0i,T O ,Ao,/Jo,r7}. 

We construct the density distribution on a three- 
dimensional grid, with spatial resolution of 0"05 = 
5.6 (t//112pc) AU, which is half the spatial sampling of the 
MICHELLE detector (0."lpix _1 ). The emission coefficient 
ju (Eq. [3]) is sampled at the grid points and numerically in- 
tegrated along each line of sight. As in the simple subtrac- 
tion in CL2] the PSF is derived from the image of HD 20893 
(Fig. [U>). We upsample the PSF to match the grid resolu- 
tion, shift it according to x„, and convolve it with the disk 
model's integrated emission. These procedures take place in 
the Fourier domain with their Fourier-equivalent operations 
via the convolution and shift theorems. We transform back to 
the image domain and bin the result to the resolution of the 
instrument (a factor of 2 in each dimension). 

3.3.2. Fitting Process 

Given the model disk, how do we characterize the range 
of model parameters allowed by the data? As we have only 
limited information (e.g. the thermal emission resolved in 
a single band), we expect that degeneracies will exist be- 
tween some parameters (such as Ao and To). Further, be- 
cause Eqs. [T]-[5] contain significant nonlinearities, we expect 
non-Gaussianity in the joint distribution of allowed parame- 
ters. Finally, this problem contains nuisance parameters (e.g. 
d) which are unimportant when considering the physical ar- 
chitecture of the disk. These characteristics, combined with 
a desire to incorporate prior knowledge, suggest a Bayesian 
approach to model fitting. The range of disk architectures 
allowed by our simple model can be described by the joint 
posterior distribution of model parameters, namely 

P (e\vj) ex P (v,Q\i) = p{v\e,r)p{Q\i). (6) 

The posterior distribution of parameters is proportional to 
the likelihood, p(T>\Q,I), times the prior, p(Q\I). Here, D 
represents the observational data and / represents our back- 
ground information. In this framework, the best-fit model 
is one whose parameters maximize the posterior distribution 
(so-called maximum a posteriori fitting). The Bayesian ap- 
proach has previous ly been used in modeling mid-IR image s 
of debris disks (e.g. lKoerner et al.lll998HWahhaj et al.ll2003l) . 
A more general comparison of Bayesian and frequentist meth- 
o ds in the contex t of common astronomical problems is given 
in lLoredol (119921) . 

Brute-force characterization of this distribution is not com- 
putationally efficient. The large number of model parameters 
ensures that the posterior has a relatively high dimensional- 
ity, and grid-based schemes for mapping the posterior require 
many evaluations of the model. Instead, we turn to Markov 
chain Monte Carlo (MCMC) methods to more efficiently ex- 
plore the parameters' joint posterior distribution. We note that 
these advanced statistical techniques stand in contrast to the 
relative simplicity of our model; our goal of measuring the 
ranges of allowed disk architectures motivates the expendi- 
ture of computational resources in the exploration of param- 
eter space rather than the complexity of the underlying disk 
model. 
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The MCMC framework uses a Markov chain, consisting 
of a series of states {0«}, to explore the parameter space 
and sample the posterior. A requirement of MCMC methods 
is an algorithm for randomly choosing states in the Markov 
chain. Given a state 0„, the Metropolis-Hastings algorithm 
provides a general-purpose method for randomly choosing the 
subsequent state 0„+i based on a candidate transition func- 
tion (commonly Gaussian). Often this algorithm is used in 
conjunction with the Gibbs sampler, which provides a pre- 
scription for choosing which of the model variables to change 
in each transition. With the Metropolis-Hasting algorithm 
choosing states, the chain will eventually reach convergence. 
A key property of converged chains is that the distribution 
of states approximates the posterior distribution. The use 
of MCMC in astrophysical data analysis has grown in re- 
cent years; a summary of this technique and a n app l icatio n 
to quantify uncertainty in orbit fitting is given by Ford (2005). 

When addressing questions about the physical properties of 
the disk, we are sometimes interested only in the range of a 
subset of parameters allowed by the model. In the Bayesian 
framework, the distribution of a parameter (or joint distribu- 
tion of a subset of parameters) is obtained through marginal- 
ization. The posterior is integrated over the variables not 
of interest. For example, to obtain the joint marginal dis- 
tribution of to and Ao, we integrate over the other variables 
e / = {rf,x*,PA i ; ) £,7,0 Ol 0i ) /«o,»7}, 

P(t ,X \VJ) = J p(Q\V,I)de'. (7) 

Having obtained an MCMC chain, this integral is trivial — 
the joint marginal distribution is simply the joint distribu- 
tion of the variables' samples, in this case the distribution of 
{to,k, Ao.„}. 

Having chosen the sampling scheme, we must define the 
likelihood and prior distributions for our model parameters. 
We calculate the likelihood by subtracting the model image 
from the observed emission presented in Fig.QJ. We compute 
xf m for the image fit using the background noise level, and 
we add the Xsed rising from the IRAS SED measurements 
(including upper limits). The likelihood is then p(V\Q,I) oc 

ex P[-(xL + XsE D )/ 2 ]- 

We have several priors whose product constitutes p(Q\I). 
We use normal priors for the Hipparcos parallax (8.92 ± 
1.05 mas; a proxy for d) and the stellar flux factor £ Q3.lt . 
We also apply a log-uniform prior to Ao, with limits of 1 nm- 
1 mm, as well as a uniform prior on 7, ranging from -4 to 4. 
The prior for the outer disk extent 0\ is taken to be uniform 
from 9q to 200". With these likelihood and prior functions, 
we ran a total of six chains using the Metropolis-Hastings al- 
gorithm with the Gibbs sampler. Each chain contained 10 4 
samples. The first 10 3 samples were discarded, as the can- 
didate transition function proposal variance for each parame- 
ter was adjusted during this "burn-in" period. We expect the 
Markov chains have converged because no gross deviations 
exist when comparing the marginal distributions of each pa- 
rameter across all chains. 

In general, the convergence rate and outcome of fitting pro- 
cesses are dependent on starting conditions. Rather than input 
hand-picked parameters to the MCMC algorithm, we first ap- 
ply a genetic algorithm (for global optimization) followed by 
a Levenberg-Marquardt least-squares fit (for local optimiza- 
tion) to find a model of the image with suitable parameters. 
During this process, the (nuisance) distance is fixed at 1 12 pc. 
Examination of the marginal distributions of each parameter 



TABLE 1 
Best-Fit Model Parameters 



parameter Model I Model II Notes 



PA 44? 5^|-| 44° 5 ±2? 5 ■■■ 

i 90°0±4.°7 90°0±4°8 ••• 

g 1.02 ±0.03 1.03 ±0.03 X7.22 x 10~ : 

7 0.07 ±0.22 -2.1+J-2 

00 0. 54 „ „ 7 0. 65_„ 09 

01 >33 7 ' > l'/S 
log 10 (r ) -«7!jj3{ -1-94^7 ... 

log 10 (Ao/lAtm) -0.65+° ; 2 i -0At%£ ■■■ 

Dependent or fixed parameters: 

h 10(<f/l 12 pc)AU 

T) 

C5 61+ 2 !J AU 73+. 2 |AU 

CJi >3700AU >210AU 

r(03o) 171+31 K 142+2; K 



NOTE. — 95% confidence intervals for marginal posterior dis- 
tributions. PA, position angle; t, inclination; £ = (R„/d) 2 , stellar 
flux parameter; 7, surface density power law index; 9q 1, disk in- 
ner/outer edges, angular units; tq vertical optical depth to absorp- 
tion at inner edge; Ao, effective grain size; ho, scale height at inner 
edge; r\, scale height power law (cf. Eq.ffJ; BJ0.1, inner/outer edges, 
spatial units; r(Cio), dust temperature at inner edge. 
il Lower-limit in Model II set by prior. 

shows that the posterior distribution in the region of the global 
maximum is smooth. As a check for other maxima, we drew 
500 randomly distributed samples of 9 over a wide distribu- 
tion, and performed a least-squares fit with each sample as the 
starting condition. All fits converged on the same maximum, 
which suggests the absence of lesser maxima in the posterior 
distribution. 

3.3.3. Results 

Using the ring model presented in 33.3.11 we fit the A/' -band 
image, and in the SED the 25 and 60 /mi points while respect- 
ing the 12 and 100 /im upper limits (Model I). Confidence in- 
tervals for the marginal posterior parameter distributions are 
given in Table [T] We are able to obtain a reasonable fit to 
the image, though the model tends to overpredict the 25 and 
100 /im fluxes (Fig. |2). With the parameters that maximize 
the posterior, the \ 2 f° r the image is 912 (with 32 2 pixels), 
while the x 2 for the IRAS SED is 3.3 (with 2 data points and 
2 upper limits). 

The marginal distributions for the position angle and incli- 
nation are approximately normal, with PA = 44.°4± l.°3 and 
i = 90.°0 ± 2.°4. The PA is marginally inconsistent with the 
direction of the inner scattered light disk, though this should 
not be taken as evidence of separate disks since we did not 
calibrate the detector orientation. We note that the range of 
allowed inclinations may be affected by our choice of disk 
scale height. 

The SED data at > 25 /im constrain the values of the outer 
radius 6\ and density power-law index 7. From the marginal 
parameter distribution, we find 6\ > 30", which is further 
than the scattered light disk has been detected (15"; iKalasI 
2005). The model also requires 7 ~ 0, which places a large 
amount of cool material in the outer disk. With Ao ~ 0.2 /mi, 
giving a m i n <C 1 /im, we can infer this to be the same pop- 
ulation giving a blue /?-[F110W] color, expected for small 
grains approaching the Rayleigh regime (iKalasI 120051) . We 
note that Rayleigh grains scatter light quasi-isotropically, in 
the sense that the first moment of the scattering phase func- 
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tion (g = (cos 8)) is zero. For an edge-on, wedge-shaped disk 
(77 = 1) of a power-law distribution of isotropic scatterers, the 
midplane surface brightness profile at projected distance b is 
oc Z? 7 " 1 . Under these assumptions, Model I predicts a bright- 
ness profile power law index of ~ — 1, which is shallower 
than the me asured indices (-2.7 to -3.7; Schn eider et al.ll2005t 
iKalasI 120051) . This suggests that the long-wavelength SED is 
not produced by a relatively flat distribution out to large radii, 
but rather by a separate population of grains. Similar popula- 
tions ha ve been invoked for o t her systems (e.g. in (3 Pic and 
AU Mic; iBackman etaill 19921; iFitzgerald etail2007l) . 

The large spatial extent of Model I and its inconsistency 
with the scattered light profile power law suggest that we 
amend the model with an additional grain population. One 
possibility would be to maintain a distribution of small grains, 
similar to Model I, to produce the A^'-band image and 25 /mi 
flux, and add larger, cooler grains that reproduce the 60 /im 
flux but contribute little to the shorter-wavelength emission. 
This would require additional parameters; modeling a single 
ring of large grains requires the ring distance and the grains' 
size and number. In the interest of computational simplicity, 
we choose not to explicitly parameterize the larger grains in 
Model II. Instead, we implicitly allow for this population by 
relaxing the SED fit, requiring only that the 60 /im flux from 
the grain population producing the spatially resolved emission 
does not exceed the IRAS measurement. The unmodeled pop- 
ulation is presumed to supplement the small grain emission to 
match the observed 60 /im flux. The larger grains are assumed 
to emit little at 25 /im. In effect, we treat the 60 /im detection 
as an upper limit. 

After running the MCMC chains with the same procedure 
as before, we find that Model II reproduces the A^'-band im- 
age and the 25 /im flux. With the parameters that maximize 
the posterior, the \ 2 for the image is 902, while the % 2 f° r 
the IRAS SED is 0.4 (with 1 data point and 3 upper limits). 
One drawback to our decision to avoid explicitly parameter- 
izing the large grains in this model is that it precludes rigor- 
ous comparison of the goodness-of-fit to Model I. 4 However, 
we expect that when combined with scattered-light data, this 
model will be favored because of its ability to include surface 
density gradients with 7 < 0. 

Model II produces less stringent requirements on 7 than 
from Model I, favoring steeper dropoff, though flat distribu- 
tions are not excluded. The distributions of the inner radii, 
in terms of both angular (#0) and physical (B3o) variables, are 
given in Figures [4^ and[4j3. 

To give a basic example of possible parameters for the un- 
modeled large grain component, we added the flux from a 
single ring of large grains at 80 AU to the best-fit Model II. 
We found that 10 29 cm 2 of Ao = 85 /im grains were able to re- 
produce to 60 /im flux without strongly affecting the 25 /im 
component from the small grains. These large grains have 
a temperature of 50 K according to Eq. [4] These large-grain 
parameters are only a single sample in a range of possibili- 
ties, and we note that they are expected to be strongly covari- 
ant without additional observations of the far-IR and sub-mm 
flux. 

In both models, the effective grain size Ao and the fiducial 
optical depth to are covariant. This is evident in their marginal 

4 This is usually done by comparing the "evidence" for each model, 
through the ratio f p(D, ei|Model l,I)dQi/ f p{T>, 6 n |Model lIJ)dBn, 
with a possible additional factor representing prior preference for one model 
over the other (e.g. Sivia 2006, Ch. 4). 
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FIG. 4. — Marginal posterior distributions of a subset of model parameters. 
Panels (a) and (b) give the probability density of location of the inner disk 
edge in angular (9q) and physical (B5o) units, respectively. The ordinates of 
these plots are matched assuming a distance of 112pc, and the differences 
evident between the distributions highlights the value of the joint estimation 
of all parameters. Panel (c) shows the joint marginal posterior distribution 
for effective grain size Ao and vertical optical depth to absorption at the inner 
edge t . 

joint distributions, shown in Figure [4};. The relaxation of the 
far-IR SED requirement allows for a broader range for these 
parameters in Model II. We stress that joint estimation of pa- 
rameters is frequently necessary in disk models. For example , 
in scattered-light modeling of AUMic, Graham et al. (2007) 
revealed a strong degeneracy between scattering asymmetry 
and the surface density power law, which was broken by the 
independent constraint from polarization measurements. 

In principle, the choice of the prior p(Q\I) can affect the 
distributions of parameters given in Table Q] A parameter's 
distribution is sensitive to the prior if it is not dominated by the 
likelihood term in Eq.|6] We explored the sensitivity to choice 
of prior by comparing the marginal distributions of each pa- 
rameter with its prior. In general, both Models I & II show 
similar sensitivity to priors. In each model, the marginal pos- 
terior distribution of £ is roughly Gaussian in shape, however 
its mean is offset from the mean of the Gaussian prior. In 
this case, both the likelihood and prior play a significant role 
in shaping the posterior distribution. The marginal posterior 
distribution of the distance d is completely dominated by the 
prior derived from the Hipparcos parallax. This is unsurpris- 
ing, since we did not use a prior for the stellar radius or lumi- 
nosity that might serve to constrain the distance. The marginal 
posterior for the outer radius 0\ largely follows the flat prior 
(out to 200"), though it has an inner cutoff set by the likeli- 
hood function. Thus the lower limits listed in Table Q] are rel- 
atively insensitive to the prior. The marginal posterior of 7 is 
well-constrained by the likelihood function in Model I, how- 
ever in Model II the prior affects the posterior. In the latter 
model, the distribution shows an upper-limit cutoff near 7 = 
and a peak near 7 = -2. However, the distribution is smooth 
down to the lower cutoff set by the prior at 7 = -4. There- 
fore the uniform prior for 7 plays a role in shaping the pos- 
terior distribution in Model II, and the likelihood function is 
only able to exclude values in the upper end of the range. For 
all other variables, the likelihood function is sharply peaked 
compared to the relatively flat priors; therefore their results 
are insensitive our choice of prior. 

4. DISCUSSION 

The spatial locations and inferred sizes of the thermally 
emitting grains are largely consisten t with the dust seen in 
scattered light. ISchneider et alJ (120051) find that the scattered- 
light disk is symmetric for 0."5 < 6 < 1 ."7, which ove rlaps 
with the symmetric mid-IR emission we measure in 33.21 The 
inner edge of the warm dust ring (6q ~ 0"5-0."7) corresponds 
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to the outer boundary of the scattered-light brightness asym- 
metry seen inward of 0"5. A comparison of the NE and 
SW F110W surface brightness profiles shows the SW ansa 
is ~ 2 times brighter than the N E at Q."5. Integra ting the disk 
light in each ansa from 8 > Q."3. ISchneider et all find the total 
NE (SW) emission to be 1.67 ± 0.57 mJy (3.14±0.57mJy), 
which corresponds to a NE-SW fractional flux difference of 
-0.31 ±0.18. In comparison to the scattered light, in the N'- 
band image we do not find an increasing brightness trend for 
the S W ansa inward of 0"6, nor do we find evidence for asym- 
metry in this region. How can the emission appear asymmet- 
ric in scattered light, but symmetric in thermal emission? One 
possibility is that the grains responsible for the scattered-light 
asymmetry reside inward of the warm ring. In this case, they 
must have properties that differ from the majority of grains 
in the inner disk, which scatter light in the symmetric com- 
ponent. If the asymmetrically distributed grains have a sig- 
nificantly higher albedo or larger average size, they may not 
be apparent in the thermal emission. Future observations of 
scattered light at different wavelengths may constrain such 
changes in grain properties with position. Another possibil- 
ity is that the asymmetry is produced by a density enhance- 
ment at sufficient distance from the star such that its grains do 
not produce significant iV'-band emission relative to the warm 
ring. 

For Ao — 0.1-0.5 fim and reasonable assumptions on grain 
composition and size distribution ( 33.3. IK we find ami n ~ 
A /6 ~ 0.02-0.1 nm. As noted in €031 in the R and F110W 
bands, submicron grains will scatter blue as their sizes ap- 
proach the Rayleigh regime, cons istent with th e blue color 
of dust scattered light inferred by Kalas (2005). This lends 
support to the contention that the scattered-light images and 
thermal image presented here are probing similar populations 
of grains. Simultaneous modeling of the scattered light and 
thermal emission, deferred for future work, has the potential 
to strongly constrain grain sizes, locations, and composition 
for the bulk of the inner disk. 

These observations are a stepping stone to understanding 
the nature of the disk and source of its structure. However, we 
lack direct evidence for the age of the system, which compli- 
cates such analy sis. Based on the stellar distance and galactic 
space motion, lKalasl(l2005h argues for an age ~30 Myr based 
on the system's possible association with the Gould Belt or re- 
cent star formation in the Taurus-Aurigae association. As we 
will show in the following calculations, this age indicates that 
the grains in the inner disk are likely not primordial and must 
be replenished. In the next subsection, we seek to estimate 
orders of magnitude for processes governing grain production 
and removal for the sizes and spatial locations in the inner- 
most disk. 

4.1. Grain Dynamics 

Assuming the radiative coupling efficiency averaged over 
the stellar spectrum Q m( i = 2, a grain density of 2 gem -3 , 
and M* = 1.8M Q , th e fiducial radiatio n pressure blow-out 
size is Gbiow ~ 3 fim (Burns e t al.lll979l) . In the absence of 
forces other than radiation and gravity, grains smaller than this 
size are removed from the system on the free-fall timescale, 
~ 10 4 yr. At first glance, this would imply that production 
of submicron grains must be extremely rapid. However, our 
conversion from Ao to the geometric size a m ; n is quite uncer- 
tain, and may be too small by a factor of ~ 30 (if the grains 
are weakly absor bing and have a size distribution steeper 
than dn oc q~ 3 5 t/a: lBackman et all l992). Another possibility, 



which maintains consistency with submicron scattered-light 
grains, is that drag forces can increase the residence time of 
the grains undergoing ejection. Finally, we note that abiow is 
not a stri ct lower limit to the steady-state size distribution. As 
noted by iBurns et all (Il979t) . very small grains (sizes much 
less than the peak wavelength of stellar radiation) couple in- 
efficiently to the radiation field, such that Q m a <C 1 and the ra- 
diation force cannot overcome gravity. For such small grains, 
composition and porosity can can play a crucial role in the 
residence ti me due to their effects o n the optical constants. 
For example. iGrigorieva et al.l d2007l) show in their Fig. 1 that 
around an A5V star, silicate grains are blown out regardless 
of size or porosity. In contrast, the smallest icy grains can 
remain bound. Compact icy grains have a blow-out size of 
3 /im, whereas «biow for grains with 80% porosity drops to 
0.7 fim. Because of their small Q rd d, the radiation pressure on 
80% porous, icy grains smaller than 0.05 /im is insufficient to 
overcome gravity. This cutoff size decreases to ~ 0.02 fim for 
more compact icy grains. We conclude that our model grain 
sizes may be compatible with the steady-state size distribution 
arising from radiation pressure, and note that this mechanism 
provides a natural separation between a population of small 
(a < 0.05 /im) grains and a population of grains larger than 

^blow- 

Collisions may play a significant role in the lifetimes 
of grains in the inner disk, as fragments from a catas- 
trophically dispersed grain can be rapidly removed via ra- 
diation pressure. For a low-eccentricity disk of single- 
size particles, the timescale between mutual collisions is 

fcoii ~ (^ T I e °) i where £1 is the orbital frequency and r^ eo 
is the geometric vertical optical depth. We find f co u ~ 

6 x 10 4 (r/60AU) 3/2 (rf "/lO" 3 )" 1 yr. Calculation of this 
timescale from our models is complicated by the unknown ab- 
sorption efficiency, as r^ eo (03o) = tq/Q^. We also note that 
the timescale for destructive collisions is likely different from 
the above f co n, and it must account for the unknown grain size 
and velocity distributions. The presence of a gas disk can 
damp relative velocities, increasing the timescale for destruc- 
tive collisions. These damping forces can also cause grains to 
settle in the mid plane, causing the d isk to be very thin in ver- 
tical extent (e.g. Gar aud et alJ20 04). Observations of the disk 
scale height may constrain the presence o f gas and its effect 
on collisional timescales. Our models in 33.31 show that the 
mid-IR image is consistent with an edge-on, vertically unre- 
solved disk. More work is needed in modeling scattered-light 
images (at higher spatial resolution) to determine if the appar- 
ent disk thickness is the result of vertical extent or inclination 
effects. Alternately, it may be possible to constrain the gas 
density through modeling of th e dust density's radial struc - 
ture, as has been done for Pic (Thebault & Aug ereaul2005l) . 
Regardless, future dynamical models of the inner disk must 
make a detailed accounting grain-grain collisions. 

Drag forces can decrease the periastra of grain orbits, fill- 
ing the inner disk with material. Assuming circular orbits, 
the Poynting-Robertson drag timescale at radius r is fpR ~ 

10 6 (a/1 /mi) (p/2gcm- 3 ) (r/60AU) 2 (ftad/2)" 1 yr. This is 
significantly longer than f co u < 10 5 yr, suggesting P-R drag 
is dynamically unimportant for the radii of warm dust emis- 
sion (though it may allow s ome mid-IR emitting dust in sys- 
tems without planets; see IWyattl 120051). The d etection of 
a potentially massive gas disk by Redfield (2007) suggests 
that gas drag may affect the grains in this system. As noted 
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above, the role of gas drag in the inner disk is difficult to es- 
timate due to our ignorance of the gas disk's density distri- 
bution and physical state, and no physical features (such as 
sharp outer edges of dust ring s or midplane settling) currently 
provide such ind ications (e.g. Takeuchi & Artvmowic j |200 It 
iBesla & Wull2007l) . 

An interesting physical result from our models are the 
warm temperatures of the effective grains at the inner rim 
of the dust annulus [T((S ) ~ 1 30-200 K]. At these temper- 
atures, the sublimation of water ice in small grains is efficient 
( f sub ^ 1 yr), suggesting no water is present in grains at this 
distance from the star. However, t his process is a strong func - 
tion of temperature (cf. Eq. 16 of lBackman & Parescdll993h . 
and therefore stellar distance (Eq.|4|l. At a distance of r = 2C3o, 
a 0.2 /im water ice grain has f su b ~ 1 Myr. This raises the pos- 
sibility that, rather than being the location of grain creation 
and outward diffusion, the inner edge of the warm dust disk 
is the destruction site for icy inspiraling grains. For this to be 
the case, drag forces must overcome radiation pressure and 
destructive collisions for a significant population of icy grains. 

The inference of the physical processes responsible for the 
disk structure (and that govern its evolution) are unclear. Di- 
rect spectroscopic evidence of the stellar age is still needed. 
Furthermore, the spatial distribution of the gas disk is an 
important direction for future observations. This system 
is attractive for the direct de tection of gas emission, simi- 
lar to observations of /3Pic dThi et al.ll2001t lOlofsson et al.1 
1200 U [Brandeke r et al.l2004h . The next modeling steps should 
combine the available data, including the resolved scattered 
light, thermal emission, and SED. More detailed calcula- 
tions can reveal the processes responsible for the disk struc- 
tur e, in a manner s i milar to that developed for AUMic 
by lStrubbe & Chiangl(l2006l) . 



4.2. Conclusions 

We have (1) spatially resolved the thermal emission from 
the warm inner disk around HD 32297 in the N' band, (2) 
found that the stellar SED is inconsistent with the tempera- 
ture and luminosity of AO; rather, we favor a cooler, less- 
luminous star, (3) found that the observed N' emission with 
the stellar PSF subtracted suggests a symmetric, optically thin 
ring model, (4) modeled the thermal emission (including IRAS 
SED data) with an annulus consisting of a single population of 
efficiently absorbing, inefficiently radiating grains, (5) deter- 
mined that a separate population producing the 60 /im emis- 
sion improves the fit to the SED and the consistency with the 
scattered light disk, and (6) identified the possibility that ice 
sublimation may play a significant role in the destruction of 
grains in the warm inner disk. 
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